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Abstract 

A covariant model based on the Bethe-Salpeter formalism is proposed 
for investigating the solar proton burning process pp De^Vf, and the 
near-threshold deuteron disintegration via electromagnetic and weak inter- 
actions. Results of numerical calculations of the energy dependence of rele- 
vant cross sections and the astrophysical low-energy cross section factor Spp 
of the proton burning process are presented. Our results confirm previous 
canonical values, and the energy dependence of the Spp factor is rather close 
to phenomenological extrapolations commonly adopted in computations of 
solar nuclear reaction rates. 



I. Introduction: The proton-proton fusion reaction pp De^v^. plays an impor- 
tant role in understanding processes occurring in the Sun and in investigating the solar 
structure. This reaction is the first step in the chain of nuclear reactions producing the 
solar energy and neutrino fluxes. Unfortunately the corresponding cross section is far 
too low for a direct laboratory study. In this context reliable theoretical calculations and 
estimates of the reaction rates are of a great importance. Nowadays, the most precise 
estimates have been obtained within the framework of low-energy potential models ■ 
The accuracy of the obtained results is believed to be of the order of few per cents or 
even better 0,3- For instance, the low-energy cross section factor Spp that determines 
the rates for the proton-proton fusion is estimated as Spp = 4.00(1 ± 0.007) x 10~^^ MeV 
b [0. While this number is considered as a canonical value, recent estimates [^], based on 
an effective relativistic field theory, claimed a considerably larger value being inconsistent 
with helioseismic data J^,!!- This has triggered a series of recalculations of Spp within 
various approaches (cf. PHTI]] and further references therein), which recover values of S.^ 



pp 



centered around the canonical one. However, there are some other indications |jT2| that 



Spp may be 4.2 x 10~^^ MeV b or even larger. 

Calculations of the pp rates require the evaluation of two main ingredients: (i) the 
weak-interaction matrix element, and (ii) the overlap of the pp and deuteron wave func- 
tions. In most relativistic effective field theories a self-consistent treatment of the deuteron 
bound state is difficult, so one usually adopts the short-range effective theory |T^ to avoid 



an explicit use of the deuteron wave function. Within such models one may achieve good 
fits of data in the zero energy limit, whereas a computation of matrix elements at finite 
values of the energy is not yet available. 

In this paper we calculate the relevant matrix elements within a fully covariant model 
of the deuteron, based on the Bethe-Salpeter (BS) formalism. Within this model a vari- 
ety of quantities, including the total and differential cross sections, energy dependences, 
angular distributions etc. are accessible. In the first part of this paper we introduce our 
formalism and consider as examples the processes of near-threshold deuteron disintegra- 
tion by electromagnetic and weak interactions. The second part extends the formalism 
and is entirely dedicated to the solar pp burning processes. Within the relativistic impulse 
approximation a fully covariant expression for the Spp factor is obtained. 

II. Deuteron disintegration near thresholds: In this section we consider processes 
of the type ID —>■ l'NiN2, i.e. the deuteron disintegration by leptons in electromagnetic 
[eD —>■ e'pn) and weak {e~^D vpp) reactions near threshold. There the amplitudes fall 
down rapidly and the cross sections are governed by the phase space of final products. 



The threshold energies are Efhr = 2.7375 MeV {eD reactions) and Ethr = 0.93139 MeV 
(weak interaction), where Ethr is the total energy of the incident lepton in the deuteron 
rest frame. We consider the electromagnetic disintegration process within the one-photon 
exchange approximation (cf. fig. la) and the weak positron-nucleon interaction (cf. fig. lb) 
within the effective Fermi model. The deuteron vertex D NN is treated in a covariant 
way within the BS formalism. The BS equation has been solved numerically |l^, with 



realistic one-boson-exchange interactions, and suitably parameterized [|T3|. The solution 
describes fairly accurate the main low-energy and static characteristics of the deuteron. 



such as the deuteron binding energy, and its quadrupole and magnetic moments |16 



The invariant cross section of the process D{1, l')NN, depicted in fig. |l], has the form 
_ (27r)^5(^) (/c; - M -PATi -PiVa) I . . ,2 dK c?p'i dp'^ 



where A(s, Mj, /i^) = [s — {Mfi + /i)^] [s — [M^ — /i)^] is a kinematical factor with s as usual 
Mandelstam variable, and fi are the deuteron and the electron (positron) masses, 
respectively. The 4-momenta of leptons are denoted as k = (£^,k), while for hadron 
momenta the notation p = {E, p) is used. The prime superscript specifies the momenta 
of the final particles. The subscript "2" denotes the spectator nucleon, i.e. the neutron 
in eD reactions and the non-interacting proton in e^D processes (see fig. |l|). The two 
nucleons in the final state are treated in the plane wave approximation, i.e. the effects of 
final state interaction are here disregarded and all calculations are performed within the 
relativistic impulse approximation. 

The square of the invariant amplitude M.id-,i'nn may be cast in the form 
\Mid~^vnn\^ = S'^{Q'^) l^f.W^'^, where S^Q"^ = -{ki - k[Y) is the Feynman propaga- 
tor of the intermediate exchanged boson (photon or boson), and the leptonic {l^u) 
and hadronic {W^u) tensors are defined by 

i,u = \y. (njMl/)(^li^lO = ^Tr [(k, + /.)r,(k; + /i')r.] , (2) 

spins 

W,. = ]-Y.{KK\J,\D){D\J+\N[N'^), (3) 

spins 

where and denote the leptonic and hadronic current operators, is the corre- 
sponding vertex function, related to the choice of interacting currents; k is used as a short 
hand notation for 7^^^^. The electron-nucleon electromagnetic vertex in the on-mass-shell 
form, F^^(Q^) = •y^Fi^Q'^) +i^^^^KF2(Q^), contains the electromagnetic form factors of 
the nucleon Fi 2, and is its anomalous magnetic moment k. In view of the small trans- 
ferred momenta and low initial energies, considered here, the second term may be safely 



disregarded, and with high accuracy also the nucleon formfactor Fi^Q"^) is replaced by 
unity. Then the electromagnetic part of the leptonic tensor is computed straightforwardly. 
For the weak vertices one may employ the Weinberg-Salam theory for weak interactions, 
however, in the considered energy range the effective Fermi model is appropriate. Thus 

3t,{x) = i)e{x)-i^ (1 - -f^)i)u{x), (4) 

X(a;)(^^=o) = cose4N{x) ^^'^'\ , (1 - Ri,)^n{x), (5) 
Gp 

and S{Q'^) = —i=, and R = 1.2681 is the ratio of axial to vector coupling constants, 
v2 

Gp = 1.16639 X 10~^ GeV~^ is the effective Fermi constant for the weak nucleon decay, 
6c = 0.24 is the Cabibbo strangeness mixing angle (AS" = means the strangeness 
conserving part of the current), iPi^n{x) are the corresponding leptonic and nucleonic field 
operators, r stands for the isospin Pauli matrices. 

With virtue of the Mandelstam technique [0] for computing matrix elements within 
the BS formalism, the deuteron tensor W^iy in eq. @ receives the form 

W^, = i (p^^ - m^)Tr [^^'ip,,p2)T^ ip[ + m,) ^^^(^1,^2) (^2 + "^2)] , (6) 

where "^^^ is the deuteron amplitude as defined in refs. [|l^,|l^,|l^ and pi and p2 = P2 
are the momenta of the deuteron constituents with a mass m, i.e. the proton (mp) 
or neutron (m„) masses. By inspecting eq. (D one observes that, due to the adopted 
mechanism in fig. |l], the spectator nucleon is on mass shell, and the first factor in eq. 
is always zero. Nevertheless, the BS amplitudes are themselves singular when one nucleon 
is on mass shell with the singularity of each amplitude being of the order {p2 — mg)^^. 
Consequently, the r.h.s. of eq. (|^) is finite. In the general case the BS amplitude \E'^"^, 
which is a 4 X 4 matrix, may be projected on a set of basis matrices and such an operation 



determines eight independent partial amplitudes |[16| , p!8[| . However not all these amplitudes 



equally contribute to the deuteron matrix elements. One may use a basis set of spin- 



angular matrices [|T8l which allows to classify the partial amplitudes corresponding to 



their contributions in a full analogy with the non-relativistic CclSG clS 

^Pi ^ ^\ ^ waves in the spectroscopical notation (the superscripts denote the p 

spin of each constituent, cf. |jl6| for details). Then the ^Sf"^ and ^Df^ waves turn out to 



be a relativistic generalization of the non-relativistic 5* and D waves in the deuteron |]TB| . 
In the considered reactions the contribution of D waves vanishes near the threshold (as 
do the other amphtudes with at least one negative p spin), leaving us with contributions 
only from the ^Si^'^ waves. A direct evaluation of traces in eqs. (0, P) results in analytical, 
covariant expressions for the cross section. For instance, the electromagnetic tensors read 



/^,iy^'^ = 16M,(27r)3$2(po,|p2|)x 

{kep[){k'pi) - /i^(pipi) - {k%)ml + {k'^p[){kePi) + 2fi^ml\ , (7) 



where ^{po, \p2\) is the covariant generahzation of the deuteron wave function ||T6[ within 
the BS formahsm, which has the following form in the deuteron rest frame 

GsiPo, IP2I) 



where Gs{po, IP2I) is the BS partial vertex corresponding to the ^Si~^ amplitude and 
m is an average value of the nucleon mass. Inserting eq. (0) into eq. (|l|) yields 

da^D = !^ $2(|p2|) dp2 Sip[' - ml) X (9) 
% J , {{keP[){k'p,) - APiP'i) - (KKH + ik'A){kePi) + 2f,'ml) ^ 



where the expression enclosed in the curly brackets together with the S function is merely 
the cross section da'^P of an electron scattering off a moving proton with the momentum 
Pi. Hence, for the electromagnetic eD disintegration processes the cross section may be 
written as 

^'""^ / dP2'^<^Mp2\)da^nk,Pi). (10) 

By processing in an analogous way with the weak disintegration diagram displayed in 
fig. lb, one obtains 

^:"^ = / dp2^^\\p2\)dal-{k,p,), (11) 

where, for convenience, the "elementary" cross section da'^{k,pi) for a hypothetical 
e~^n up reaction has been introduced by 

da::^ik,p,) = {K + 1)) \ikePi){kM + {keP\){kM] 

+ (R^ - l))m,mn{kX) - [{kep\){k'^Pi) - {k,p^){k'^p\)]] (12) 

In computing eq. ([llD the angular integration on p2 should be performed in only one-half 
of the hemisphere due to the two identical particles in the final state. 

In fig. results of the numerical calculation of the total cross section as a function of 
the initial energy of the electron (positron) in the laboratory frame are displayed. Above 
the electromagnetic disintegration threshold there is probably no way to experimentally 
separate the weak processes from the electromagnetic ones. Below the electromagnetic 



and above the weak thresholds, the cross section of positron weak disintegration of the 
deuteron seems far too low to be measured under laboratory conditions. So one may 
conclude that at present an experimental investigation of the matrix element, relevant 
for the solar pp fusion, is still impossible also when using the cross channels. (This is in 
contrast to other reactions, cf. [^].) Only reliable theoretical calculations may be applied 
to estimate the nuclear reaction rates and to interpret the corresponding experimental 
neutrino data. 

III. Solar proton burning process: Let us now consider the process of the two- 
proton fusion P1P2 ve^D at low relative energies. While this reaction looks rather 
similar to the processes in the previous section, here appears a peculiarity. Due to the 
relation 2mp > + fi (with nip = 938.27231 MeV, = 1875.61339 MeV, deuteron 
binding energy Ed = 2.22455 MeV, and fi = 0.510999 MeV) there is no energy threshold 
for the pp fusion and, in principle, the fusion process may occur at arbitrarily small relative 
energies of the two protons. In this case, in the cross section (|I]) two factors will play the 



major role: (i) the two-proton flux (~ y^A(s, m^, m^)) approaches rather rapidly zero and 
the computation of the cross section becomes problematic, (ii) a strong Coloumb repulsion 
between protons makes the amplitude 7W sharply decreasing and compensates singularities 
appearing from the vanishing flux. In order to avoid uncertainties connected to these 
circumstances one usually separates the contribution of Coulomb repulsion by factorizing 



the Coulom b par t of the two-proton final state, which at low relative energies may be taken 
Q^v - 1 



as ippp = ~ y'^'^V 6 ^'^ ■> where 77 = a/f is the Sommerfeld penetration parameter. 



V denotes the relative velocity of the protons, and Tpp = ■nipv'^/A is the corresponding 

relative kinetic energy; the exponent is known as the Gamov penetration factor p|. Then 

one introduces the notion of the Spp factor, which is determined only by the nuclear 

S (T ) 

part of the interaction and is connected with the cross section by cr = — pp_^-2-Kr) ^ 
Then the rate of pp fusion reactions can be written in the form |I|J^ {av)pp = 1.3005 x 
10~^^ {2 /Tq)^^^ f Spp{Tpp)e~'^ cm^s~^; here Tq is the temperature in units of 10^ K, r ~ 
15 — 40 dominates the temperature dependence of the reaction rate, / is a screening factor 
calculated first by Salpeter [^, and Spp is in MeV b. The Spp factor is needed at the 



most probable interaction energy Tj(°), which is in a range of 5 keV ■ ■ ■ 30 keV. However, 
in most analyzes in the literature, the value of Spp is quoted at vanishing energy, obtained 
, e.g. within the low-energy effective theory PJ^,p!3[1. Hence, for the sake of extrapolation 
to the finite value of Tj^^ a good estimate of the associated derivatives is also required. 

We now apply the same formalism as in the previous section, i.e. the relativistic impulse 
approximation, to estimate the low-energy behavior of Spp. The effect of initial state 



interaction is accounted for only in the Coulomb part of the interaction by adopting 
the corresponding part of the above two-proton wave function tppp. A treatment of the 
nuclear part of the proton wave function as plane waves means that in the matrix element 
all transitions are allowed. Nevertheless, at low energies only the lowest partial waves 
contribute to the proton wave function and one may expect that the total cross section 
will be determined by the axial part of the current operator and correspondingly by 
the super-allowed 0+ 0+ Gamov- Teller transition. 

In the limit Tpp —>■ 0, apart from the phase space volume, the value of the deuteron 
wave function '?/'(Po; |p|) at vanishing values of its arguments becomes important. As 
can be seen from eq. (H), at |p| 0, 2Ep — * 2m the wave function ip{po, |p|) displays 
a sharp maximum. Our solution of the BS equation |T^ has been obtained by solving 



a system of coupled integral equations with Gaussian meshes with a minimum value of 
IpI = 0.310848493733 MeV, so that an extrapolation to |p| = needs to be employed. 
This procedure causes some numerical uncertainty which we estimate to be about 10% in 
the final results. 

By observing that in the center of mass of colliding protons 2^X{s, nip, rrip) = 
4E1E2 V = sv and adopting the Gamow penetration factor, eq. ([1|) leads to 

= §^ I^^I^^ \^pp-e..nf Sikl, - (13) 

where, due to the antisymmetrization of the initial two-proton state, the matrix element 
J^pp^e+uD is now determined by two Feynman diagrams, as depicted in fig. ^, with a 
relative minus sign. Correspondingly, the hadron tensor W^y consists of a direct and an 
interference term, 



W^^f ^ = (pI - ^p)Tr [^{pi,pn)r^{p2 + mp)r,^(pi,p„)(p2 + nip 

+ {pj - mJ)Tr [^'(j9i,p„)r^(pi - mp)Tf,^{p2,Pn)iPi - "^p)] , (14) 

W^^u'^ = (pI - ml){pl - ml)TT [^{p2,Pn)r,^{puPn)r, + ^{p,,Pn)r,^{p2,Pn)r.] , (15) 

where is the vertex corresponding to the effective current operator (^). The direct 
part of the hadron tensor has exactly the same form as for the reactions e^D vpp (cf. 
eq. (^)), resulting in an "elementary" cross section like that in eq. ([I2|) . The interference 
term is different and much more cumbersome, and for the sake of brevity we do not display 
it here. 

In fig. H the energy dependence of Spp computed by eq. (|13D is displayed. The Spp 
factor, starting from Tpp -|- up to solar energies T^^, displays a rather weak energy 
dependence. At very low kinetic energies the value of Spp remains fairly constant, having 



a value of 3.960 x 10^^^ MeV b in agreement with previous results 0,^. At larger energies, 
however, the Spp factor increases rather rapidly. The two curves in fig. § correspond to 
two different methods of numerical calculations. In one case the deuteron wave function 
is replaced by its value at zero transferred momenta. To some extent this corresponds to 
calculations within low-energy effective potential models. The full line reflects results of 
calculations without such an approximation. Only at relatively high values of the kinetic 
energy, Tpp > 20 KeV, the effect of the deuteron wave function becomes sizeable. 

Within low-energy effective theories the associated derivatives of the Spp factor play an 
important role. The best estimates of the first derivative seem to be those quoted in ref. 
^, S'pp = 4.48 X 10-24 b, and the logarithmic one S'^^/S = 11.2 MeV'^ In figs. | and | 
we display our estimates for the energy dependence of these derivatives. The derivatives 
are computed at each value of Tpp by interpolating Spp with a quadratic interpolation 
formula. From these pictures it is seen that within our approach both the derivative and 
the logarithmic derivative are below the estimates in [Q] by ~ 15 — 20%. However, at 
energies around T^^^ our results for the Spp factor and the results of an interpolation with 
parameters from are rather close to each other. 

IV. Summary: We present a covariant model, based on the Bethe-Salpeter formalism, 
to investigate the energy dependence of the matrix elements relevant to the solar proton- 
proton burning process. Explicit covariant expressions for the cross section and the low- 
energy cross section factor Spp are presented. Our numerical estimates of Spp and its 
derivatives confirm previous results. 
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FIG. 1. Feynman diagrams for electromagnetic (left panel) and weak (right panel) dcuteron 
disintegration. The weak coupling constant and the propagator of the boson has been 
replaced by the effective Fermi constant Gp for the 4-point contact interaction. 
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FIG. 2. The threshold behavior of the total cross section in electromagnetic and weak dis- 
integrations of the deuteron. The arrows depict the corresponding thresholds. The dashed line 
displays the integrated eD enp cross section when the outgoing electron polar angle O is 
larger than 15°. 
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FIG. 3. Feynman diagrams contributing to the processes pp e^uD. 
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FIG. 4. The Spp factor for the reaction pp e^uD as a function of the relative kinetic energy 
of incident protons. The full line is the result of a calculation by eq. (^), the dotted line is the 
result of a computation with the deuteron wave function taken as constant, ^{po, |p|) = ^'(0, 0). 
The inset shows Snn vs. Tr,,, in linear scales. 
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FIG. 5. The first derivative of the low-energy nuclear reaction cross section factor Spp. No- 
tation as in fig. ^ 
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FIG. 6. The logarithmic derivative of the low-energy nuclear reaction cross section factor 
Spp. Notation as in fig. ^ 
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